Relativistic coupled-cluster single-double method applied to alkali-metal atoms 
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extensions to include a dominant class of triple excitations) led to accurate predictions for energies, 
transition amplitudes, hyperfine constants, and other properties of monovalent atoms. Further 
progress in high-precision atomic structure calculations for heavy atoms calls for improvement of 
the linearized coupled-cluster methodology. In the present work, equations for the single and double 
excitation coefficients of the Dirac-Fock wave function, including all non-linear coupled-cluster terms 
that contribute at the single-double level are worked out. Contributions of the non-linear terms to 
energies, electric-dipole matrix elements, and hyperfine constants of low-lying states in alkali-metal 
atoms from Li to Cs are evaluated and the results are compared with other calculations and with 
precise experiments. 
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I. INTRODUCTION 

A relativistic version of the coupled-cluster single- 
double (CCSD) approximation for monovalent atoms is 
developed. In this approximation, single and double ex- 
citations of the (frozen-core) Dirac-Fock wave function 
for an atom with one valence electron are included to all 
orders in perturbation theory. The relativistic CCSD is 
an extension of an earlier relativistic single-double (SD) 
approximation, in which all nonlinear coupled-cluster 
terms were omitted. The relativistic SD approximation 
provides a method closely related to the configuration- 
interaction method for including classes of correlation 
corrections to Dirac-Fock wave functions to all-orders in 
perturbation theory. When modified to include the dom- 
inant triple excitations, the SD method led to accurate 
predictions for energies, transition matrix elements, life- 
times, hyperfine constants, and polarizabilities in alkali- 
metal atoms (see, for example, Refs. [E S H S S Hi)- 
Owing to recent improvements in the precision of atomic 
experiments, it is important to go beyond the relativistic 
SD approximation and include the non-linear coupled- 
cluster terms. 

The coupled-cluster method is an all-order extension 
of many-body perturbation theory introduced into nu- 
clear physics by Coester and Kiimmel ^, [1] and adapted 
to atomic and molecular physics by Cizek and Paldus 
[§, [13, [HI- A comprehensive review of the coupled- 
cluster method and it s ap plications in quantum chem- 
istry is given in Ref. [l^. Relativistic CCSD calcula- 



tions for monovalent atoms have been carried out previ- 
ously by Ynnerman et al. Il3l for transition energies in 
Li-like U, by Eliav et al. [14| for ionization and excita- 
tion energies of alkali-metal atoms, by Lim et al. flS^ for 
polarizabilities of alkali-metal atoms from Li to element 
119, by Martensson-Pendrill [l6l | fo r charge radii of fran- 
cium isotopes, by Landau et al. l3| for electron affinities 
in alkali-metal atoms Na to element 119, by Chaudhuri 
et al. [l^ for ionization and excitation energies in Rb 
and Cs, by Porsev and Derevianko 1^ for properties 
of Na, by Derevianko and Porsev [20 1 for properties of 
Cs, and by Sahoo et al. plj for parity nonconservation 
in Cs and Ba+. Among the nonrelativistic applications 
of the CCSD method to alkali-metal atoms, we mention 
those of Lindgren [S^l who calculated energies and hyper- 
fine constants of 2s and 2p states in Li and Salomonson 
and Ynnerman "23^ who evaluated energies and hyperfine 
constants of 3s and 3p states in Na. 

The relativistic SD method for alkali-metal atoms was 
introduced in Ref. [1^, where it was used to obtain en- 
ergy levels, fine-structure intervals, hyperfine constants, 
and electric-dipole matrix elements in Li and Be"*" . Later, 
the SD method was applied successfully to study proper- 
ties of heavier alkali-metal atoms and monovalent ions 
[E S, i, i, i, l- Despite the success of the SD all- 
order method (and its extension to include important 
triples) in predicting many properties of monovalent sys- 
tems, there are various instances where it fails to produce 
accurate results. For example, the magnetic-dipole hy- 
perfine constant for the Id^ /2 state of Fr calculated using 
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the SD method significantly disagrees with the experi- 
mental value [25j]. The ndc^i'^ hyperfine constants in Cs 
provide other such examples [2a]- This issue is not lim- 
ited to the hyperfine constants; correlation corrections 
to electric-dipole Qdj — 6pji matrix elements in Rb are so 
large that they change the sign of the matrix elements [3] ■ 
The large differences between third-order and all-order 
values noted in Ref. [1] lead to the conclusion that the 
omitted higher-order contributions are significant. Cases 
where correlation corrections are extremely large lead to 
poor accuracy for the SD method, which is generally of 
high-precision; in such cases, a more accurate treatment 
of correlations is mandatory. 

Another motivation for further development of the rel- 
ativistic SD method is the study of parity nonconserva- 
tion (PNC) in heavy atoms. One goal of PNC studies is 
to provide atomic-physics tests of the standard model 
of the electroweak interactions through determination 
of the weak charge Qw] another is to extract nuclear 
anapole moments from PNC measurements. A precise 
calculation of the PNC amplitude (together with an un- 
certainty estimate of the calculation) is necessary to de- 
rive a value of the weak charge Qw from experimental 
measurements. The accuracy of the most advanced ex- 
perimental study of PNC (which was carried out in Cs) is 
0.3% [27], while the accuracy of the corresponding theo- 
retical calculations is about 0.5% (see HI] and references 
therein). The difference between the value of Qw ex- 
tracted from experiment and the value inferred from the 
standard model stands at la More accurate exper- 
iments, coupled with improvements in the calculation of 
PNC amplitude will lead to a more accurate determina- 
tion of Qw and a more stringent test of the standard 
model. 

Recently, a rather complete treatment of a class of 
triple excitations that included the valence electron and 
two core electrons was carried out in Refs. [l^ [13] for 
Na and Cs, respectively. Contributions of quadratic va- 
lence non-linear terms were also calculated and found to 
be relatively large. As a result, further development of 
the all-order method must include a complete treatment 
of the non-linear terms at the SD level. 

In this work, we include all valence and core non-linear 
coupled-cluster corrections to the SD equations, leading 
to a complete set of CCSD equations, and we study their 
effects on various atomic properties of alkali-metal atoms. 
In particular, we calculate ionization energies and hyper- 
fine constants for the ground ns states and the npi/2, 
np3/2 excited states of Li, Na, K, Rb, and Cs. Reduced 
electric-dipole matrix elements for the ns — and 
ns — np^/2 transitions are also calculated. We give a de- 
tailed breakdown of various non-linear contributions in 
order to identify the most important effects and to mea- 
sure the influence of core non-linear terms. Comparison 
of our calculations with other theoretical calculations and 
with experimental results is presented. 



II. COUPLED-CLUSTER METHOD 

In coupled-cluster theory, the wave function of an atom 
with one valence electron in a state v is written as 

l^-,) =exp(5)l$,), (1) 

where \^v) is the lowest-order atomic wave function for 
atomic state w, which taken to be a frozen core Dirac- 
Fock (DF) wave function, and where the wave operator 
(the operator that maps the DF wave function onto the 
exact wave function) is expressed in exponential form 
exp (S) . The exponential function in Eq. ([1]) can be ex- 
panded to give 

]*,) = (l + S'+i^^^---)!*.). (2) 

The cluster operator S is expressed as a sum of n-particle 
excitations Sn of the lowest-order wave function 

S = Si + S2 + --- . (3) 

As the number of excitations increases, the complexity 
of the wave function increases. The computational com- 
plexity rises dramatically beyond the double excitation 
term 5*2. Retaining only single and double excitations, 
Eq. ^ may be written 

1^,) = (1 + 5i + 52 + Isf + S1S2 + Isf 
2 

+ lsl + lsfS2 + ^^St + ---)\<f.). (4) 

The one-electron excitation Si may be either an exci- 
tation of a core electron or an excitation of the valence 
electron. Correspondingly, the single core and valence 
excitations are given by 

Sic — ^ Pma^^^a^ 
ma 

Sly = ^ ^ PmvO^lri^v- (5) 

Similarly, for double core and valence excitations 

S2C — 2 ^ ^ PmnabO,m^n^b^aj 
mnah 

S2V = PninvbalnO-liO'bav (6) 

mnb 

The expansion coefficients pma and pmv are referred to 
later as single core and valence excitation coefficients, 
while Pmnab and Pmnva are referred to as double core and 
valence excitation coefficients, respectively. In Eqs. (O 
and ([5]), a\ and are creation and annihilation oper- 
ators for an electron state i. Here and in subsequent 
formulas, the indices from the beginning of the alphabet 
a, 5, ... designate core states, indices from the middle of 
the alphabet to, n, ... designate excited states, the index v 
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labels the valence state, and indices i, j, k and I designate 
arbitrary states. 

In the SD method, only terms linear in the excita- 
tion coefficients are retained, and all remaining terms are 
omitted. Substituting Eqs. (O and ([6]) into Eq. Q and 
retaining the terms linear in the excitation coefficients 
yields the single-double (SD) all-order wave function 



a^abaa 



mnab 



ra^v 7nnb 



(7) 



To derive the equations for the excitation coefficients (see 
Ref. [23| for details), the SD all-order wave function ([7]) 
is substituted into the many-body Schrodinger equation 



(8) 



where the Hamiltonian H is the relativistic no-pair 
Hamiltonian [2£?|, which can be written in second- 
quantized form a,s H = Hq + V, where 



Ho = y^e^ajoi] 



V 



i 

ijkl 



(9) 



2 ^9ijki[a\a\aiak\ 



(10) 



In the no-pair Hamiltonian, contributions from negative- 
energy (positron) states are omitted. Products of oper- 
ators enclosed in brackets, such as [aja]a/afe], designate 
normal products with respect to a closed core. The quan- 
tity U in Eq. (fTO|) is the model potential used in the Dirac 
equation defining single-particle orbitals. In this work, 
Uij is taken to be frozen-core Dirac-Fock potential 



(VbF)ij - y^(ffmja - 9iaaj)- 



Such a choice of the potential significantly simplifies cal- 
culation since the second term in Eq. (jTU]) disappears. 
The quantities Qijki are two-body Coulomb matrix ele- 
ments: 

g.,u = jd\jd\'4{v)4>]{v')y^^M^)i>i{v') 

and the quantity ti in Eq. ^ is the eigenvalue of the 
Dirac equation. The third term in (jlOp is a c-number 
and provides an additive constant to the energy of the 
atom. 

The all-order SD wave function given by Eq. ([7]) in- 
cludes only terms that are linear in the excitation coef- 
ficients. In the present work, we take into account all 



non-linear terms that arise from the single and double 
excitations. Out of all possible non-linear terms, only six 
terms, \Sl, S1S2, \Sl, \Sl, \SlS2, and ^Sf contribute 
to the single-double equations. Explicitly, the non-linear 
terms contributing to the core single-double equations 

are 2*^ic' ^icS2c^ q^ic^ 2*^20' '2^'ic^^c'> and 24*^ic and 
the non-linear terms contributing to the valence single- 
double equations are SicSiy, {81^820, SicS2v} , \8i^8iv, 
82c82v, {8iv8ic82c, ^8i^82v}, and ^8f^Siy. 

First, we consider the contributions from the non- 
linear core terms. The first three non-linear core terms 



^•Sic = PrcPsdalalaaOc, (11) 

rscd 

SicS2c ^ 2 ^ ptePrscdalalalaeadac: (12) 

rstcde 

7:8fc = ^ X! PtePrcPsdaWsalaeOdac (13) 



contribute to equations for both single and double exci- 
tation coefficients, while the last three terms 

^4 = ■^'S'L = o X! PrscdPtuefalalala'lafaeadac,{14:) 



rstucdef 



T5 - Ti^lc^- 



2c 



T ^ PtePiLf Prscdalalalal^afO, 



(15) 



rstucdef 

^6 = -^'S'^c = ^ PrcPsdPtePufala\a\aiafaeadac 



rstucdef 



(16) 



contribute to the equation for the double excitation co- 
efficients only. 

The right-hand side of the single-double equations is 
obtained by operating on the non-linear terms above with 
the two-particle interaction operator 



G = ]^'Y9i3ki[a\a]aiak] 

ijkl 



(17) 



according to Eq. pH)) . 

To derive the equation for the core single-excitation 
coefficients, we extract those terms in GTu {k = 1, 2, 3) 
that are proportional to aJ^OaoJ. To derive the equa- 
tion for the core double-excitation coefhcients, we extract 
the terms in GT^ {k = 1 • ■ • 6) that are proportional to 
^al^al^Qhaaal. In all cases, we drop terms correspond- 
ing to disconnected diagrams. Below, we use the no- 
tations GT^ and GTj^ to designate contribution of the 
corresponding terms to the single or double excitation 
equations, respectively. For clarity, we give the contribu- 
tions from all terms separately. 

The equation for the core single-excitation coefficients 
becomes 



SD + GT^ + GT^ + GT^ 



3 ' 



(18) 
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where Ci is the one-body DF energy for the state i, SD is 
the contribution from the linear coupled-cluster terms 
given in [l^j and the contributions of the non-linear 
terms are 

GTi — gmdrsPraPsd — QcdasPmcPsd i 

(19) 

drs cds 

GTi = - Y.~9cdsrPrsdaP^c~Y.~9^asrPs^cdPra 

cdrs cdrs 

H~ ^ ^ QcdrsPrmcaPsdi (^0) 
cdrs 

QcdsrPmcPrdPsa ■ 

(21) 

We used the notation gmnab = Qmnab - gmnba and 
Pmnab = Pmnab - Pmnba in the abovc formulas. 

The equation for the core double-excitation coefficients 

is 

(Ca + 66 ^ Cm ^ ^n) Pmnab — SD 

+ GT^' + GTa'' + GTi + GT^ + GT^ + GT^, (22) 

where 

GT^ — ^ ^ fJninrsPraPsb ^ ^ ^ ScdabPmcPnd 



^ ^ dmdarPrbPnd 



a h 
m ^ n 



(23) 



^ ^ QcdrbPndPrmca ^ ^ QcdarPrdPmnch 
cdr cdr 

^ ^ 9cdraprhPnmcd H" ^ ^ ^ 



ncr s Prh P smca 



cdr 



^ ^ QncrsPscPrarab ^ ^ QncrsPmcPsrah 
crs crs 



a ^ b 
m ^ n 



(24) 



GT3'' 



^ ^ Qcdar PndPmcPrb ^ ^ QracrsPncPraPsb 
. cdr crs 



a ^ b 
m n 



(25) 



GT4'' 



^ ^ Qcdtui^ptubdpninac ~(~ PmucdPntha) 
cdtu 

a ^ b 
m ^ n 



GT^ — ^ 9cdtu{PtaPub + 



PmcPndPtuab ) 



cdtv 



.cdtu cdtu 



QcdtuPtbpncPmuad) + 



cdtu 



a ^ b 



GT^ = E ScdtuPtaPubPmcPnd- 



(27) 



(28) 



All non-linear contributions to double-excitation coeffi- 
cients are symmetrized to preserve the property Pmnab — 

Pmnba • 

Only one non-linear term, GTi = ^GS^c, contributes 
to the equation for the core correlation energy: 



bmnPinaPnb •, 



(29) 



abmn 



where (J-Ef is the core correlation energy obtained with 
linearized SD wave function ([7]): 



-\Y.3ab 



'ibmn Pmnab • 
mnab 

We note that the summation over each index, for ex- 
ample i, involves summing over the principal quantum 
number n^, the relativistic angular momentum quantum 
number Kj, and the magnetic quantum number m^. The 
sum over the magnetic quantum numbers is carried out 
analytically and the final formulas are given in Appendix 

The equations for the valence excitation coefficients 
Pmv and Pmnvb are identical to the core equations given 
by Eqs. (|18ll28p with replacement of index a by index v 
and addition of the valence correlation energy 6Ey into 
the parenthesis on the left-hand side of both equations, 
i.e. 

(e, - em + SE,)pmv = SD + [GT^ + GTi + GT^)^^^ , 

(30) 

(e„ -I- £6 — Cm — e„ + SE^) Pmnvb = SD + (^GTf 

+GTi + GTi + GTi + GTi + GTi) . (31) 
The valence correlation energy SE^ is given by 

5Ey = SEf-Y, Qcdvtptdpvc \ / Qvdtuptvpud 
cdt dtu 

^ ^ 9cdutPvcPutvd ~ ^ ^ QcdutPtv Puvcd 
cdtu cdtu 

+ Qcdtupvtvcpud — dcdutPtdpuvpvc- (32) 



cdtu 



cdtu 



The term SE^^ represents the expression for the valence 
correlation energy without the non- linear terms [2^ : 

(26) 5i?.'°=E 

QvavmPnia^ / QabvYnPinvab^ / QvbmnPmnvb- 
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The non-linear contributions to the valence correlation 
energy arise from the SicSiy, {SiyS2c, SicS2v}, and 
^Sf^Siv terms. 

We solve the SD equations using a finite basis set. Each 
orbital wave function is represented as a linear combina- 
tion of the B-splines. We consider a radial grid of 250 
points within a sphere of radius 100 a.u. We include 35 
out of 40 basis orbitals for each angular momentum and 
include all partial waves with I < 6 in our calculations. 
A detailed description of the B-spline method is given 
in Ref. [lOl- We treat the non-linear terms on the same 
footing with the linear terms, i.e. all the linear and non- 
linear terms are iterated together. First, the equations 
for the single core and double core excitation coefficients 
are iterated until the core correlation energy given by 
Eq. (P^j) converges to relative accuracy e = 10~^. Then, 
the valence equations are iterated until the valence corre- 
lation energy given by Eq. (|32p converges to the relative 
accuracy e. Atomic properties can be evaluated once the 
values of the excitation coefficients are known, as briefly 
described below. 

Matrix elements of a one-body operator Z = 
Zija\aj are determined using the formula 



(33) 



Substituting the expression for the wave function from 
Eq. ([7]) in the above equation and simplifying, one finds 



v/(i + iv„)(i + iv^)' 



(34) 



where Zc, 



Ny, and iVy, are linear or quadratic 



functions of the single and double excitation coefficients 
written out in Refs. [ll.[23|. 

In general, the non-linear terms coming from expand- 
ing the exponent in the CC wave function also contribute 
to the expressions to matrix elements. Even at the CCSD 
truncation level, one encounters an infinite number of 
such contributions. A rigorous method of partial sum- 
mation (dressing) of the resulting series was devised in 
Ref. (3l| . The method is built upon expanding a prod- 
uct of cluster amplitudes into a sum of n-body insertions. 
Although in the present paper we do not include these 
direct non-linear contributions to matrix elements, calcu- 
lations [20] show that dressing may contribute as much 
as a few 0.1% to hyperfine constants in Cs. 



TABLE I; Contributions of the non-linear (NL) terms to re- 
moval energies for Li, Na, K, Rb, and Cs. A detailed descrip- 
tion of all contributions is given in text. SD designates the 
correlation corrections to the energies calculated using the SD 
method. All results are in cm~^. 



T \ 




2pi/2 


2P3/2 


SD 


406.0 


352.1 


352.0 


Core NL terms 


0.9 


0.1 


0.1 


S2cS2v 


-5.1 


-4.2 


-4.2 


SlcSlv, {SlvS2c, SlcS2v} 


-2.6 


-0.9 


-0.9 


(Jtner valence iNL bu terms 


0.0 


0.0 


0.0 


Total 


399.2 


347.1 


347.0 


iNa 


■3S1/2 


3pi/2 


3P3/2 


SD 


1483.3 


462.0 


459.8 


Core NL terms 


0.5 


2.6 


2.5 


S2cS2v 


-44.3 


-15.3 


-15.3 


SlcSlv, {SlvS2c, SlcS2v} 


-23.8 


-10.2 


-10.1 


(Jtner valence INL ou terms 


0.0 


0.0 


0.0 


Total 


1415.6 


439.1 


436.9 


Iv 


4Sl/2 


4pi/2 


4p3/2 


SD 


2869.7 


1114.3 


1100.5 


Core NL terms 


26.8 


11.3 


11.2 


S2cS2v 


-142.5 


-59.1 


-58.5 


SlcSlv, {SlvS2c, SlcS2v} 


-62.9 


-34.4 


-34.3 


(Jtner valence INL bu terms 


0.1 


0.1 


0.1 


Total 


2691.2 


1032.2 


1019.1 


riu 


OSl/2 


OPl/2 


5p3/2 


SD 


3423.2 


1301.1 


1236.3 


Core NL terms 


31.7 


13.5 


13.1 


S2cS2v 


-185.4 


-76.7 


-73.2 


OlcOlv, XUlv02c, JlcJ2v f 


n 

- lUO.U 


4-7 Q 


-4:U.O 


Other valence NL SD terms 


0.3 


0.1 


0.1 


Total 


3164.8 


1190.1 


1130.0 


Cs 


6Sl/2 


6pi/2 


6p3/2 


SD 


3881.5 


1618.7 


1442.3 


Core NL terms 


44.3 


18.1 


16.8 


S2cS2v 


-224.4 


-107.3 


-96.3 


SlcSlv, {SlvS2c, SlcS2v} 


-162.4 


-74.6 


-68.3 


Other valence NL SD terms 


0.9 


0.4 


0.4 


Total 


3539.9 


1455.3 


1294.9 



conducted five separate calculations for each atom. Each 
subsequent calculation includes all terms in the previous 
calculation together with additional terms, the effect of 
which is being determined. For clarity, we describe each 
of the calculations below. 



III. RESULTS AND DISCUSSION 
A. Energies 

Table |T] shows a detailed breakdown of contributions 
from non-linear terms to the removal energies of the 
alkali-metal atoms Li, Na, K, Rb, and Cs. To illus- 
trate the relative importance of the various terms, we 



1. Linearized SD calculation, with all non- linear terms 
omitted. The results of this calculation are listed 
in the rows labeled "SD" . 

2. All non-linear terms are included in the core equa- 
tions only, all non-linear valence terms are omitted. 
The differences of those values and the SD results 
are listed in rows labeled "Core NL terms" . 

3. All non-linear terms are included in the core equa- 
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tions and the quadratic term S2cS2v is included in 
the valence equations. The differences with the cal- 
culation (2) give the contributions from the non- 
linear valence term S2cS2v and are listed in the rows 
labeled accordingly. 

4. All non-linear terms are included in the core equa- 
tions, and all remaining quadratic terms are in- 
cluded in the valence equations. The differences 
between those values and the results of the calcu- 
lation (3) give the contributions of the SicSiv and 
{SivS2c, 5'ic'5'2ti} quadratic valence terms. 

5. Final calculation: all non-linear terms are included 
in the core and valence equations. The results of 
this calculation are listed in the rows labeled "To- 
tal". The differences between those values and 
the results of the calculation (4) give the contri- 
butions of the cubic and quartic non-linear terms 
that are listed in rows labeled "Other valence NL 
SD terms". 

In all of the cases considered here, the addition of 
the non-linear terms results in a decrease in the corre- 
lation contributions to the removal energies. From Table 
m we see that while the contribution of the core non- 
linear terms to the removal energies is negligible com- 
pared to contribution of the valence non-linear terms for 
the ground states of Na, it becomes significant (over 10% 
of the total NL contribution and about 1% of the total 
correlation energy) for all states of K, Rb, and Cs con- 
sidered here. Furthermore, the contribution from core 
NL terms is opposite in sign to that from the valence 
non-linear contribution. 

The S2cS2v term gives the dominant nonlinear contri- 
bution, as expected. More than half of the contribution 
from non-linear terms to the removal energy arises from 
this term. However, contributions from other quadratic 
terms, SicSiy and {81^820, SicS2v}, are also signifi- 
cant. Finally, the contributions from terms ^Sf^Siy, 
{SiySicS2c, ^Sf^S2v}, and ^Sf^Siy, i.e. terms that are 
cubic or quartic in the excitation cocfhcicnts are negli- 
gible. For Li and Na, contributions from these terms 
are essentially zero; for K, Rb, and Cs, the contribu- 
tions are less than 0.5% of the total from all non-linear 
terms. The breakdown of contributions from the valence 
NL terms is essentially identical for all states of Na, K, 
Rb, and Cs under consideration. The total contribution 
of the NL terms to the correlation energies of the lowest 
three states of Na, K, Rb, and Cs is remarkably large; it 
is about 1.5% of the linear SD correlation energy for Li, 
5% for Na, and 9% for Cs. 

In Table |TT1 we compare of our results for the corre- 
lation with the linearized all-order Sp(pT) calculations 
of Ref. [2], CCSD calculations of Ref. 0], the CCSD(T) 
calculations of Ref. [l^ and with the experimental ener- 
gies • The values in the rows labeled "SD" are the 
sum of lowest-order (DF) energies, the SD contributions 
given in Table U and the extrapolated contributions of 



TABLE II: Comparison of the all-order removal energies with 
the theoretical results obtained by the coupled-cluster method 
^14, .18] and experiment. All the values are in cm~^. 



Li 


2Sl/2 


2pi/2 


2p3/2 


SD 


43494 


28586 


28585 


CCSD 


43488 


28581 


28580 


Expt. [32] 


43487 


28584 


28583 


Na 


3Si/2 


3Pl/2 


3p3/2 


SD 


41444 


24495 


24476 


CCSD 


41376 


24472 


24453 


SDf2] 


41447 


24494 


24477 


CCSD [14] 


41352 


24465 


24447 


Expt. [32] 


41449 


24493 


24476 


K 


4Si/2 


4pi/2 


4P3/2 


SD 


35258 


22219 


22066 


CCSD 


35080 


22044 


21984 


SD[2] 


34962 


22023 


21966 


CCSDfl4] 


35028 


22016 


21957 


Expt. [32] 


35010 


22025 


21967 


Rb 


5Sl/2 


5pi/2 


5p3/2 


SD 


34021 


21241 


20994 


CCSD 


33762 


21130 


20888 


SD[2] 


33649 


21111 


20875 


CCSD[14J 


33721 


21117 


20877 


CCSD[18] 


33603 


21080 


20831 


Expt. [32] 


33691 


21112 


20874 


Cs 


6Si/2 


6pi/2 


6p3/2 


SD 


31871 


20421 


19842 


CCSD 


31529 


20258 


19695 


SD[2] 


31262 


20204 


19652 


ccsDyj 


31443 


20217 


19669 


CCSD [18] 


31250 


20137 


19574 


Expt. [32] 


31407 


20228 


19674 



the higher partial waves (Eoxtrap); the values in the rows 
labeled "CCSD" are the sums of the DF, SD, non-linear 
contributions, and Eoxtrap- The CCSD values should 
agree with the calculation of Ref. [3l within the nu- 
merical uncertainties of the calculations. Significant dif- 
ferences between all-order SD results and CCSD results 
were noted earlier in Ref. 0, indicating that NL terms 
may be large. In this work, we find a good agreement be- 
tween our CCSD values and the results of Ref. [14| . The 
remaining discrepancies can be explained by the differ- 
ences in some numerical details of the calculations. In our 
calculations, we truncated the number of partial waves 
included in all sums over excited states at /max = 6. Con- 
tributions from partial waves with I > 6 are extrapolated 
in second order (see Q for details of the extrapolation 
procedure). The resulting correction increase for heavier 
alkali-metal atoms and is about 0.1% of the total energy 
of 6s state in Cs. The differences between our values and 
those of Ref. [ij are of the same order of magnitude as 
the Eoxtrap contributions. 
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The SD approximation omits contributions to ener- 
gies of third order in perturbation theory that arise from 
triples excitations. These missing third-order terms were 
included in the calculation of Ref. The size of this 
terms is approximately given by the differences between 
values hstcd in rows "SD" and "SD"2]". The calcula- 
tion of Ref. ^] omits these terms entirely. Interestingly, 
these terms are of nearly the same magnitude as the NL 
terms and of the same sign. As a result, both calculation 
of Ref. 2], that omitted NL terms, and of Ref. [U], that 
omitted missing third-order terms, were in quite good 
agreement with experiment. It was shown in [Tol l2Clj|. 
that actual iterated triple contribution significantly dif- 
fers from the third-order values. For this reason, and to 
make a clear comparison with previous CCSD results, we 
omit all triples contributions in this calculation. For Rb 
and Cs, we also list values from the CCSD(T) calcula- 
tions of Ref. ^ in the rows labeled "CCSD[li|" that 
include contributions from triple excitations, but omit 
the odd-parity channels. 



B. Reduced electric-dipole matrix elements 

Table IIIII gives the detailed breakdown of the contri- 
butions of the non-linear terms to the reduced electric- 
dipole matrix elements for Li, Na, K, Rb, and Cs. This 
table is structured in exactly the same way as Table [H 
the only exception is that we added lowest-order DF 
value to the linearized SD values in rows labeled "SD" . 
The breakdown of the non-linear terms is identical to 
the breakdown for correlation energies. The major con- 
tribution comes from the term S2cS2v as in the case of 
the removal energies. Also, there is almost no contribu- 
tion from terms which are cubic or quartic in the exci- 
tation coefficients. The core non-linear terms and the 
cubic and quartic valence non-linear terms decrease the 
El reduced matrix elements, while the quadratic valence 
non-linear terms increase these values. The only excep- 
tion to the statement above is the contribution of the 
core non-linear terms for Na which behaves in the op- 
posite manner. Since the contribution of the quadratic 
valence non-linear terms outweigh the contribution of all 
the other terms, the ultimate effect is an increase in val- 
ues of the reduced electric-dipole matrix elements. As in 
the case of the correlation energies, the contribution of 
the NL terms is rather large, especially for heavy alka- 
lies where it reaches 1.5%. The comparison of our results 
with the experiment [H, |33| is also given in Table IIIII 
Addition of non-linear terms results in a deterioration of 
the agreement with experimental values as expected from 
owing to significant cancellation between the va- 
lence NL and valence triple terms. A very interesting 
finding of this work is relatively large contribution of the 
core non-linear terms for heavy alkalies. While the core 
non-linear term is entirely negligible for Na, it is found 
to be 0.2% for Cs. 



TABLE III: Contributions of the non-linear terms to the 
nsi/2 — npi/2 and nsi/2 — np3/2 reduced electric-dipole ma- 
trix elements for Li, Na, K, Rb, and Cs. The final values are 
compared with experimental results. All values are given in 
atomic units (eao, where ao is the Bohr radius). 



Li 


2Si/2 - 2pi/2 


2Sl/2 - 2p3/2 


SD 


3.31654 


4.69033 


Core NL terms 


-0.00007 


-0.00008 


S2cS2v 


0.00063 


0.00088 


OlcOlv, \ Jlv02c, JlcJ2v f 


0.00020 


0.00028 




n nnnnn 


n nnnnn 

u.uuuuu 


Total 


3.31730 


4.69141 


Expt. [33] 


3.317(4) 


4.689(5) 


Na 


3Si/2 — 3pi/2 


3Si/2 — 3p3/2 


SD 


3.53099 


4.99314 


Core NL terms 


0.00005 


0.UUU06 




0.00487 


0.00690 


cc rcc cci 

OlcOlv, \ JlvJ2c, JlcJ2v f 


0.00211 


0.00297 




n nnnnn 

u.uuuuu 


n nnnnn 

u.uuuuu 


Total 


3.53802 


5.00307 


Expt. [33] 


3.5246(23) 


4.9838(34) 


K 


4Sl/2 - 4pi/2 


4Sl/2 — 4p3/2 


SD 


4.09820 


5.79392 


Core NL terms 


-0.00474 


-0.00669 


S2cS2v 


0.02261 


0.03198 


OlcOlv, \Jlv02c, JlcJ2v f 


0.00873 


0.01231 




n nnnn9 

-u.uuuuu 


n nnnn'^ 


Total 


4.12478 


5.83149 


Expt. [33] 


4.102(5) 


5.800(8) 


Rb 


5Sl/2 - 5pi/2 


5Sl/2 — 5p3/2 


SD 


4.22005 


5.95527 


Core NL terms 


-U. 00649 


-0.00913 


S2cS2v 


0.03230 


0.04564 


CC rcc cci 

OlcOlv, \iJlv'J2c-, >JlciJ2v f 


0.01533 


0.02156 


Other valence NL SD terms 


n nnno/i 

-U.UUUU4 


-U.UUUUD 


Total 


4.26115 


6.01328 


Expt. [33J 


4.231(3) 


5.977(4) 


Cs 


6Sl/2 - 6pi/2 


6Sl/2 - 6p3/2 


SD 


4.48157 


6.30391 


Core NL terms 


-0.01057 


-0.01482 


S2cS2v 


0.04585 


0.06437 


SlcSlv, {SlvS2c, SlcS2v} 


0.02762 


0.03865 


Other valence NL SD terms 


-0.00015 


-0.00021 


Total 


4.54432 


6.39190 


Expt. [34] 


4.4890(65) 


6.3238(73) 



C. Hyperfine constants 

We calculate the magnetic-dipole hyperfine constants 
A for the nsi/2 ground states and the ex- 
cited states of Li, Na, K, Rb, and Cs. The nuclear mag- 
netization density is described by a Fermi distribution 
with half-density radius c and 90%-10% falloff thickness 
t = 2.3 fm. Table Hvl lists values of the parameter c(fm) 
used for magnetization distribution and the gyromag- 
netic ratio gi for each of the alkali-metal atoms. The 
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TABLE IV: Values of the Fermi half-density c(fm) parameter 
used for magnetization distribution and gyromagnetic ratio 
gi (in units of the nuclear magneton) which were used to 
calculate hyperfine constants. 



Atom 


c 


gi 


^Li 


1.7995 


2.17093 


23Na 


2.8853 


1.4784 


39k 


3.6108 


0.260993 


»5Rb 


4.8708 


0.54136 


"3CS 


5.6748 


0.737886 



S'l term was found to be dominant for removal energies 
and electric-dipole matrix elements with other quadratic 
terms being also significant. In the case of the hyperfine 
constants, the contributions of the other quadratic terms 
exceeded that of the 5*1 term for most cases. The con- 
tribution of the cubic and quartic terms was found to be 
negligible in all cases considered in the present work. In- 
clusion of non-linear terms in the single-double all-order 
method is a significant step toward further development 
of high-precision methodologies for the calculation of the 
atomic properties. 



nuclear magnetic moments were taken from Ref. [35j and 
weighted averages were considered when more than one 
value was present. 

Contributions of the various non-linear terms to the 
hyperfine constants are given in Table jVl The compari- 
son of our results with experiment [H, [33, [H, [H, Ho, IH!, 
IH,!!^ is also given. The structure of Table IVl is identical 
to that of Table Uni The "SD" values are the sum of the 
lowest-order DF values and correlation correction calcu- 
lated in the absence of the NL terms. Addition of the 
non-linear terms resulted in a decrease in the values of 
the hyperfine constants in comparison to the linearized 
SD values. The NL core terms and the cubic and quar- 
tic valence non-linear terms contribute with a positive 
sign to the hyperfine constants, while the quadratic va- 
lence non-linear terms contribute with negative sign for 
K, Rb, and Cs. For Na and Li, all core and valence non- 
linear terms contribute with the same sign. The most 
significant difference between the NL contributions to the 
hyperfine constants and to the removal energies and El 
matrix elements is that the contribution of the S2cS2v 
terms is smaller than the contribution of the other three 
quadratic valence non-linear terms S'ic5'i^, {81^820, and 
SicS2v} for the ground states of Na, K, Rb, and Cs. The 
breakdown is also diiferent for the ns, npi/2, and 71.^3/2 
hyperfine constants while it is very uniform for all ener- 
gies and El matrix elements considered here. The con- 
tributions of the core NL terms are particularly large for 
the hyperfine constants, almost 1%(!) of the total values 
for Cs. 



IV. CONCLUSION 

We have extended the relativistic SD method to in- 
clude all non-linear terms at the SD level. The effect 
of the non-linear terms on the removal energies, hyper- 
fine constants, and electric-dipole matrix elements of the 
alkali-metal atoms from Li to Cs was systematically in- 
vestigated. In particular, five different calculations were 
carried out to establish the importance of the various 
contributions for each alkali-metal atom. The effect of 
the core non-linear terms was found to be not negligible 
for heavier alkalies, reaching nearly 1% of the total values 
of the Cs hyperfine constants. Among other terms, the 
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APPENDIX A: ANGULAR REDUCTION 

1. Designations and definitions 

The Coulomb interaction gmnab is decomposed into the 
product of a term Jk (mnab) 

Jkimmab) = x 

f im k ia\f in k jb \ 

\ -mm -q rria J \ -nrin q mb J ' ^ 

depending on only the angular momentum quantum 
numbers ji and rrii of the four states (to, n, a, 5), and 
a term Xk {mnab) depending only on the principal quan- 
tum numbers rii and angular quantum numbers ki of the 
states: 

g„inab ^^Jk{mnab)Xk{mnab). (A2) 

k 

Here 

Xk{mnab) = (-1)'' (k„ ||C''|| k,,) {kji «&) Rk{mnab). 

(A3) 

The quantities Rk{mnab) are (relativistic) Slater inte- 
grals and (kto IIC**^ II 'Sa) is a reduced matrix elements 
of a normalized spherical harmonic. The quantities 
Zk{m,nab) are given by 

Zkimnab) = Xk{mnab)+y^[k] \ \xk'{mnba), 

^ { Jn Jb K j 
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TABLE V: Contributions of non-linear terms to the magnetic-dipole hyperfine constants A(MHz) of Li, Na, K, Rb, and Cs and 
comparison with experimental results. The experimental values are from Ref. 36] unless noted otherwise. 



Li 


2Si/2 


2pi/2 


2P3/2 


ou 




40.1 ( D 


-z.zy 1 


Core NL terms 


-0.025 


-0.005 


0.006 


S2cS2v 


-1.439 


-0.212 


0.167 


SlcSlv^ {SlvS2c^ SlcS2vy 


-0.183 


-0.002 


-0.011 


Other valence NL SD terms 


0.000 


0.000 


0.000 


Total 


393.585 


44.957 


-2.129 


Expt. 


401.75 


46.17(35)*''^ 


-3.07(13)''=' 


Na 




3pi/2 


3P3/2 


ou 


ooo.ZoD 


yo.uou 


1 Q OPI/l 

lo.o04 


Core NL terms 


-0.663 


-0.242 


-0.057 




-2.645 


-1.022 


-0.193 


SlcSlv^ {SlvS2c^ 5lc52ii} 


-4.878 


-0.696 


-0.174 


Other valence NL SD terms 


0.001 


0.000 


0.001 


Total 


880.101 


93.090 


18.431 


Expt. 


885.8 


94.44(13)('') 


18.534(15)'''' 


K 


^-^1/2 


^Pl/2 


4/^3/2 


SD 


237.159 


28.696 


6.214 


Core NL terms 


0.457 


0.138 


0.032 


52c52t; 


-2.783 


-0.798 


-0.123 


5ic5it;, {5lti52c, S'lcS2v'\ 


-3.208 


-0.489 


-0.145 


Other valence NL SD terms 


0.007 


0.002 


0.001 


Total 


231.632 


27.549 


5.979 


Expt. 


230.85 


28.85(30) 


6.09(4) 


Rb 


.5 '> 1 / o 
^-^1/2 


.5 1 / r, 
'-''51/2 


.571o /r, 

^yj,/2 


SD 


1051.554 


125.624 


25.560 


Core NL terms 


5.319 


0.710 


0.155 


52c52t; 


-14.198 


-4.091 


-0.577 


5ic5it;, {5iii52c, SicS2v^ 


-22.636 


-3.058 


-0.816 


Other valence NL SD terms 


0.047 


0.007 


0.002 


Total 


1020.086 


119.192 


24.324 


Expt. 


1011.9 


120.7(1) 


25.029(16) 


Cs 


6Si/2 


6pi/2 


6p3/2 


SD 


2439.053 


311.138 


51.900 


Core NL terms 


20.455 


2.604 


0.475 


52c52u 


-39.441 


-12.871 


-1.137 


SlcSlv, {5lt,52c, 5ic52d} 


-80.300 


-11.259 


-2.489 


Other valence NL SD terms 


0.321 


0.043 


0.015 


Total 


2340.088 


289.655 


48.764 


Expt. 


2298.2 


291.89(8)'^' 


50.275(3)'»' 



'"' Reference 
'''' Reference 
'■'=' Reference 
'''' Reference 
Reference 
'■'^' Reference 
'^' Reference 



where [k] = 2fc + 1. Double excitation coefficients have 
the same angular structure as Coulomb matrix elements: 

k 



and the quantities Skimnab) are defined in the same 
way as Zkijnnab). The angular reductions for the single- 
excitation coefficients are defined as follows: 
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Pmv ~ *^Km.«^u *^"^m"T't; *^(^^)' 



(A4) 



where n is the relativistic angular momentum quantum 
number defined as 



« = T(i + l/2)forj = /±l/2. 



2. Angular decomposition of terms contributing to 
the equation for core single excitation coefficients 



GTl = ^j'p\:Zo{mdrs)S{ra)S{sd) 

drs V ^-^"J 

~ H\i TAzo{cdas)S{mc)S{sd) 

[ja\ 



cds 



GTl = - ^ ^ ^?m'i° Zi{cdsr)Si{rsda)S{mc) 



- E 



c4rsl W^'J 
(^—\'^c-\-s-a- 



cdrsl 



- Zi{cdsr)Si[smcd)S{ra) 



'[3 A 



+ 4 / YT-T6j^j^ZQ{cdrs)So{r'mca)S{sd) 

\ja\ 



cdrs 



GTi = -Y,\lT^Mcdsr)S{mc)S{rd)S{sa) 

cdrs V 



GTi = ^iJ^}ll^Zi{cdrb)Si{rmca)S{nd) 



cdr 



~ E \/ MSjj^Zo{cdar)Si{mncb)S{rd) 
cdr V ^-^aJ 





I 


' \ 








d 


m J 




c n j 



cdrjk 

X Xj{cdra)Sk{nmcd)S(rb) 

— — -^TTj Zi(ncrs)Si{smca)S{rb) 

crs 

+ W Y^5j^j^ZQ{ncrs)Si{mrab)S{sc) 



X Xj{ncrs)Sk{srab)S{mc) + 





/ .n 




0| 


1^ n 


s b / 


[ TO 


raj 



a ^ b 



GTi = ^Xi{cdar)S{nd)S{mc)S{rb) 

cdr 

- ^Xi{mcrs)S{nc)S{ra)S{sb) + 



a ^ b 
m ^ n 



3. Angular decomposition of the terms 
contributing to the equation for core double 
excitation coefficients 



GTt = ^Xi{mnrs)S{ra)S{sb) 

rs 

+ ^Xi{cdab)S{mc)S{nd) 



cd 



^Zi{mdar)S{rb)S{nd) + 

_ dr 



a f-^ b 
m <r^ n 



GT 



cdtu kik^kzk 



m t C 
k2 k I 





ks 




f k2 








u 




1 TO 


a 





n b u 



Xfc^ (cdtu)Sk2 {tuab)Sk3 {mncd) 
Zi {cdtu) Si {mtac)Si [undb) 



{-^Y+d+t+u 

irTPl ^]bO<: {cdtu)Sk {tubd)Si (mnac) 



cdluk 



(-1) 



c+d-\-t+u 



- Sj^j^ Zk (cdtu) Sk {mucd) Si {ntba) 



+ 



cdtuk 

a b 
m n 
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I 




( k I 






d 




y a c 


TO J 









I 






U b j 


1 TO 


t 





GTi = ^ (_;L)n+m+a+&j; 

X Xj{cdtu)Sk{mncd)S{ta)S{ub) 

X Xj(cdtu)Sk{tuab)S{mc)S{nd) 

1 / ■p^(Sj^j^Zo(cdut)S';(TOnao?)S'(t&)5(wc) 

+ \ j^Sj^j^^ZQ{cdtu)Si{muab)S{tc)S{nd) 
cdtu V ^-^'"J 



cdt-u 



[I] 



a ^ b 
m ^ n 



GTi = YXi{cdtu)S{ta)S{ub)S{mc)Sind) 



SE, 



NL 




Zo{cdvt)S{td)S{vc) 



J f^Zo{vdtu)S{tv)S{ud) 

dtu V ^-^"J 

, Zk{cdut)Sk{utvd)S{vc) 

cdutk ^^J^-^^J 

, ,! . , Zk {cdut)Sk (uvcd) S (tv) 

cdutk ^''^ 

\ T~iZo{cdtu)So(vtvc)S(ud) 



cdtu 



cdtu 



y X T4zo(cdut)S{td)S(uv)S(i 



4. Angular decomposition of terms contributing to 
the core and valence energies 



6E^^ = - V [ja] [jb] Zo {abmn)S {ma)S{nb) 
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